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Abstract 

We present a strategy to control the evolution of a quantum system. The novel aspect of this 
protocol is the use of a single auxiliary subsystem. Two applications are given, one which allows for 
state preservation and another which controls the degree of entanglement of a given initial state. 
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I. INTRODUCTION 



Protocols for control and manipulation of quantum systems are essential for the develop- 
ment of quantum information theory [1]. Advances of this area can provide proper tools to 
avoid decoherence and to conduct quantum evolution to desirable results. Some examples 
of those strategies are: Quantum Zeno Effect (QZE) [2-4], Super Zeno Effect [5], strong 
continuous coupling [6, 7], Bang-Bang control [8, 9], etc. 

In the present contribution we show a new possibility to control the evolution of a quan- 
tum system through successive interactions with a single auxiliary subsystem. These inter- 
actions are described as unitary evolutions (in a finite time period), that can stop or inhibit 
the evolution of the system of interest. 

To make our results concrete, we present the protocol based on two examples. In the 
first one, we study a system of interest composed by two coupled qubits (Sb and S c ) sharing 
one excitation. A third two level system (S a ) is the auxiliary subsystem responsible for the 
control. Interactions between Sb and S a are inserted in the evolution of the system S& — S c 
and control its dynamics. The quantity of such interactions, as well as their duration are 
parameters that allows for several forms of control. We performed an analytical calculation 
for the state vector of the global system after N interactions with S a . The Hamiltonian form 
of this system permits the mapping of the global evolution (S a — Sb — S c ) in the real euclidian 
subspace, suggesting a geometrical interpretation for this dynamics. We also pointed out 
the differences between the present dynamics and the QZE. 

In the second example we show how to control the entanglement dynamics presented 
in Ref. [10], through interactions with a single auxiliary subsystem. In Ref. [10], two 
initially entangled atoms undergo different time evolutions. One of them interacts with 
an electromagnetic mode in a cavity and the other one evolves freely. The dynamics is 
nondissipative and the entanglement oscillates. The introduction of an auxiliary subsystem, 
that interacts with the atom in the cavity, allows for the control of entanglement dynamics. 
An empirical implementation for this process may be realized with the experimental setup 
used in Ref. [11], where a two level atom interacts with two electromagnetic modes preserved 
in the same microwave cavity. In such empirical implementation, the two level atom and 
one of the modes in the cavity (M n ) compose the system of interest, and the second mode 
(Mb) acts as an auxiliary subsystem. 
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II. TWO QUBITS DYNAMICS 



Let us consider the system of interest composed by two coupled qubits (Sb — S c ) and 
another qubit (S a ) as an auxiliary system. The hamiltonian that governs the interaction 
between Sb and S c is given by: 

H bc = e a \l a )(l a \ + e 6 |l 6 )(l 6 | + e c |l c )(l c | + I a ® hG bc (a b _a c + + o\o c _), (1) 

where a + = |1)(0|, <j_ = |0)(1|, Gb c is the coupling coefficient, I a is the identity matrix on 
the subsystem S a . The coefficients e a , e& and e c are the eigenvalues of the free hamiltonian. 

The goal is to control the dynamics in subsystem S b — S c through interactions between 
the auxiliary qubit S a and Sb- The hamiltonian for these auxiliary interactions is 



H ah = e a \l a )(l a \ + e b \l b )(l b \ + e c |l c )(l c | + hG ab (a a _a b + + o\o h _) ® I c , (2) 

where I c is the identity matrix on S c and G a b is the coupling coefficient between S a and Si,. 

Suppose that S a — Sb — S c share one excitation and e a = = e c . Since the operators H bc 
and H a b preserve the excitation number, we can write in one excitation subspace the time 
evolution operators as 



1 
cos 9 —i sin 9 
— ism 9 cos 9 



(3) 



and 
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(4) 



where 9 = G bc tbc, <P = G a bt a b, t bc (tab) is the interaction time between S b and S c (S a and S b ). 

The control of S b — S c dynamics is induced by unitary operators (XJ ab ) inserted iV times 
in the free evolution of S b — S c . The number of interventions and the durations of each one 
are the parameters that specify the control. The general expression for the state vector of 
the global system submitted to this control is given by: 



\M = (u a6 (0)u fec (^)J |^(o)) , 



(5) 



where the time evolution of S& — S c was divided by N interactions with the auxiliary sub- 
system. 

In the appendix we calculate the vector state \iPn) for a general initial state. Now, let us 
consider the initial state \ip(0)) = |0 o ) |0b, l c ) which goes through a quantum transition when 
submitted to the time evolution U& c (|) 1^(0)) = |0 o ) 1 1&, C ). We show the inhibition of this 
transition through a sequence of unitary interactions with the single auxiliary subsystem. 

As it is shown in the appendix, the global evolution of this system can be mapped on R 3 , 
therefore we may represent a sequence of N interactions with the auxiliary subsystem as: 



R 3 (0) Ri 



7T 



2N 



N 



f(0) 



ac ( 1 — cos Np) + b sin Np 
be (1 — cos Nip) — a sin Np 
(1 - c 2 ) cosiV> + c 2 



(6) 



where ^ 2rc and r(0) 



. Taking the limit N — > oo we have 



lim r/v = f(0), 

N^oo 



(7) 



as a — > 0, 6 — > and c — > 1. 

To give a geometrical interpretation of this effect consider a vector r in the euclidian 
subspace. In the present dynamics, the rotations Ri (—9) around the axis Ox are clockwise. 
Therefore, when f has a positive y component, the rotations Ri {—&) will reduce the z 
projection of the vector, but when r has a negative y component the rotation Ri (—9) will 
do the opposite, tending to compensate the previous decrease. The rotations R3 (0) move 
the vector through the subspaces where the y component is positive and where it is negative. 
Therefore, when we study the dynamics of r inserted by rotations R 3 (0) we notice that the 
decreasing of the z projection, induced by Ri (—9) when y > 0, is compensated by the 
increasing of the same projection, also induced by Ri (—9), but when y < 0. Choosing the 
angle of the rotations R3 (0) and the number of interventions N it is possible to preserve the 
projection z of r*or even freeze the dynamics of r. In Fig.l we explicitate another geometric 
point of view for the effect: the net effect of Ri (—0) and R3 (0) concerns a rotation around 
the vector n, when the relation 9/<f> decreases, n gets closer to e 3 and the state will be always 
closer to the initial one. 
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FIG. 1: The curves at the unitary sphere give the evolution of the terminal point of the vector f 
(represented with its initial point at the origin) starting at r = e% and subjected to N rotations 
concerning R3 ((f)) R4 (— ir/2N) = R(n, (p). The line segments explicitate the direction of each 
rotation axis h. (a) 4> = n/16 and three values for N: N = 10, N = 20 and N = 40. Higher N 
corresponds to n closer to £3 and to curves closer to the initial point, (b) N = 20 and three values 
for 4>: cf> = it/32, (f) = n/16 and (f) = tt/8. Higher 4> corresponds to n closer to £3 and to curves 
closer to the initial point. 



For the initial vector |-0(O)) = |0 a )|0 6 , l c ), which corresponds to (r(0)) T 
survival probability is 



(0,0,1), the 



Poor = |^v-r(0)| 2 . (8) 

The structure of hamiltonians H bc and H ab conserves the state vector in the subspace 
{|l a , 0ft, C ) , — % |0 a , 1ft, C ) , |0 a , 0ft, l c )}, allowing us to assume the presented geometrical in- 
terpretations. However, the control of quantum state through interactions with a single 
auxiliary system is not restricted to systems with such symmetry. 

From a broader point of view, this effect occurs when the accumulation of interactions 
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with the same auxiliary subsystem changes the signal and reduces the absolute value of the 
quantum transition rate. In Fig. 2 we show the function of the system S a — S& — S c . 



The oscillations of induce oscillations on the behavior of Pooi(^) (increase-decrease). 
Therefore, for appropriate values of the parameters (N and interaction time with the aux- 
iliary subsystem) we may have the increase of the function Pooi(^); at some time intervals, 
compensating the decrease in other time intervals, inducing (in average) preservation of the 
initial state, as it is shown in Fig.2. 
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FIG. 2: Probability Pooi(*) an d quantum transition rate ^ffM*) with = ^, # = ^ 

This inhibition of quantum transition induced by the increase on N is similar to the 
discrete QZE, but is structured differently: the QZE as presented in Ref. [12], has a system 
of interest interacting with N auxiliary (probe) subsystems. After each interaction the 
complete information about the occurrence of the quantum transition is available on the 
probe. This fact implies the cancelation of the transition rate after each interaction [13]. 
Therefore, we may characterize the interactions between system and probe as a measurement 
process (pre-measurement). The net effect of iV interactions between a single auxiliary 
system and the system of interest can not be characterize as a measurement process at all. 
Consequently for the later dynamics, the quantum transition rate is not necessarily null 
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after each interaction, as it is shown in Fig. 2. 

Another difference between these two effects is that in the dynamics presented here a 
transition like |l a )|0&, C ) — > |0 a )|0&, l c ) is intermediated by Sb, i-e., if only one auxiliary 
subsystem interacts iV times with the system of interest, the excitation present in S a may 
return to the subsystem S c inducing on P oi larger values than the survival probability 
observed in the QZE. A comparison between these two probabilities is shown in Fig. 3. 
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FIG. 3: Comparison between the probability of non-transition related to QZE Z(n) = cos (^y) 



and P(n) = Pqqi, with 



§ and 6 = ^ 



III. CONTROL OF THE ENTANGLEMENT DYNAMICS 



Let us consider a system composed by two space-separated atoms, one of them is isolated 
(atom B) and the other (atom A) is coupled with an electromagnetic mode (Mi), as studied 
in Ref. [10]. If there is an initial entanglement between atoms B and A (or between atom B 
and mode Mi) it changes over time, even if the initially entangled systems are not coupled. 
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This entanglement dynamics takes place because of the coupling between atom A and mode 
Ml. In this section we proposed a control of this entanglement dynamics through successive 
interactions with a single auxiliary subsystem. 
Suppose the initial state is 

\1>(0)) = \g a )(a\l 1 ,g b )+/3\0 1 ,e b )), (9) 

where the mode M x and the atom B are entangled. The coefficients a and (3 give the 
intensity of the initial entanglement. 

The coupling between atom A and Mi can be described by the Jaynes-Cummings model. 
After a time evolution the vector state of the system has the form 

\^{t)) = a(cas(gt)\g a )\l x ,g b ) -isin(^)|e a )|0 1 ,^))+/3|0 1 ,e 6 ), (10) 

where g is the coupling coefficient between atom A and Mi. We consider the atomic tran- 
sition frequency of atoms A and B resonant with the frequency of Mi. Notice that when 
t — | we have the entanglement swap, the entanglement initially present in the subsystem 
Mi-atom B is completely transferred to atom A-atom B subsystem. 

To quantify and study the entanglement dynamics we write the concurrence [14] between 
Mi and atom B (details of this calculation are in Ref. [10]). 

C Ml At) = 2\af3cos(gt)\. (11) 

The concurrence oscillates assuming null values when gt — where k is an odd number. 

It is possible to control the entanglement dynamics through N interactions with a single 
auxiliary subsystem. In the presentation of such control we consider that Mi is in a cavity 
that supports two ortogonal modes (Mi and M 2 ). An empirical realization of this system 
(in a microwave cavity) is reported in Ref. [11]. 

The modes have different frequencies (uj\— uo 2 = 5). The difference between atomic energy 
levels may be controlled by Stark effect. The detuning between the modes allows for the 
atom A, when coupled to Mi (M2), not to interact with M2 (Mi). Therefore it is possible 
to control the coupling time between atom A and the modes Mi and M 2 . 

A sequence of interactions between the atom A and M 2 , inserted into the time evolution 
shown in (9) and (10), is responsible for the control on the dynamics of entanglement between 
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atom B and Ml. The evolution of the global system is composed by N steps, each one in 
two stages. At the first stage the atom A interacts with Mi and at the second stage the 
atom A interacts with M 2 (the second stage is responsible for the control). 
In the first stage the time evolution is governed by the hamiltonian 



Hi = hwia{ai + hw 2 ala 2 + hui\e a )(e a \ + hwi\e b )(e b \ + hg(a^ai + a a a\), 

= H[ + hwi\e b ){e b \, (12) 

and in the second stage by the hamiltonian 

H 2 = hjJia\ai + ftu) 2 a\a 2 + ^2|e a )(e a | + fkui\e b )(e b \ + hg(ia^a 2 - ia^ a\), 

= H' 2 + fiwi\e b ){e b \, (13) 

where a£ = \ek)(gk\, cr^ = \gk)(ek\ (k = a,b). H[ and H' 2 act only on the subsystem 
composed by atom A, Mi and M 2 . Notice that the coupling coefficient in the second stage 
(in H 2 ) is imaginary, this is due to the ortogonal mode's polarization. 
The unitary time evolution operators for the first and second stages are: 



(14) 



Written on the basis {|0i, g a , li), |0i, e a , Oi), \li,g a , Oi)} the operators e lH i*/ R and e lH 2 t / h 
assume the form: 



e isti o 
cos(gti) —ism(gti) 
— ism(gti) cos(gti) 

cos(gt 2 ) — sin(gt 2 ) 
sin(^t 2 ) cos(gt 2 ) 



(15) 











3 — i<5*2 



(16) 



The time evolution is given by 



-iH 2 t 2 /h e -iH 

and the initial vector state of the global system is 



(17) 
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\^(0)) = \g a )(a\l 1 ,g b ) + P\0 1 ,e b ))\0 2 ), (18) 

where M 2 (auxiliary subsystem) is prepared in the vacuum state. 

After N steps the interaction time between the atom A and Mi is T = Nt 1 . Let us 
consider T = Nti = ^ (as it is in Fig.4). When there is no participation of the auxil- 
iary subsystem in the dynamics the concurrence between the mode Mi and the atom B is 
null, because the entanglement is completely transferred to the subsystem atomA-atomi?. 
Therefore, with the intervention of the auxiliary subsystem the entanglement dynamics is 
inhibited, i.e. the enhance on the number of interactions with the auxiliary subsystem allows 
for the preservation of the concurrence initial value, even when the total time for the inter- 
action of the system (atom A, atom B and Mi) is T = Nt\ = ^ (time of the entanglement 
swap), as it is shown in Fig.4. 

To summarize, we have presented a new strategy to control the evolution of a quantum 
system. This strategy requires only unitary interactions between the system of interest and 
a single auxiliary subsystem. We discuss two examples for the strategy application. In the 
first one the auxiliary subsystem controls excitation exchange between two qubits. In the 
second example, a single auxiliary subsystem is used to control the entanglement of a system 
composed by an isolated atom and a Jaynes-Cummings atom. 
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APPENDIX: CALCULATION OF \%j> N ) 

In this section we calculate explicitly the vector state \i^n)- Let us focus on 
the action of the unitary matrix XJ ab (<f))XJ bc {9) on the vector |£) written in the basis 
{|l a ,0 6 ,0 c ),-i|0 a ,l 6 ,0 c ),|0 a ,0 6 ,l c )}: 





6 




£i cos <j) — (£2 cos 9 + £3 sin 0) cos (j> 


tu(0) ij bc (o) 10 = u a6 (0)u bc (fl) 






£1 sin + (£ 2 cos + £3 sin 0) cos 
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£ 3 cos 9 — ^2 sin 
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FIG. 4: Concurrence between mode Mi and atom B with 5 = 8 x 10 5 s , g = 1.5 x 10 4 s 1 and 

gt 2 = § ■ 



This action may be mapped as a rotation on the real euclidian subspace, by choosing con- 
veniently the rotation matrix and assuming that £i, £ 2 and £ 3 (components of are real. 
Defining 



and the vector 



Ri (if) = R(ei,vj) 



R 3 (if) = R(e 3 ,f) 



1 

cos <f — sin if 
sin<^ cos ip 

cos ip — sin if 
sin ip cos f 
1 



6 

6 



(A.2) 



(A.3) 



(A.4) 



we notice that the action of the matrix R x (—9) and R 3 (<p), over r produces the same effect 
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on its components as the action of Uft c (#) and U a fc(0) on the components of i.e., 



R 3 (0) Ri(-0)r 



£1 cos — (£2 cos + £3 sin #) sin 
£1 sin + (£ 2 cos 9 + £3 sin 0) cos 
£ 3 cos # — £2 sin 



(A.5) 



The ortogonal matrix R 3 (0) Ri (—9) may be written as 

R 3 (0) Ri (—9) = R (n, if) = exp {^ph ■ J j , 



(A.6) 



where n ■ J is the generator of rotations around the axis defined by the unitary vector 
n = aii + be 2 + ce 3 , ip is the angle of rotation around h, J = e±Ji + e 2 J2 + e 3 J 3 where Ji, 
J2 and J 3 are the generators of rotation around the axis Ox, Oy and Oz, respectively. 



(A.7) 
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In order to calculate the action of (fj a b((f))~(j t, c (9) j over the state vector |^ ), we use 
the mapping of the operator U ab (0)U bc (0) on the ortogonal matrix R 3 (0) R x (—9) and the 
identity: 



[R 3 (0) Ri (—9)] N = exp(ym-j) " = exp (N(ph ■ j) 



(A.8) 



The calculation of [R 3 (0) Ri (—9)] N is reduced now to finding the axis n and the angle 



ip. 



a. Finding h and ip 
Explicitly, 



iph ■ J = ip 



-c b 
c —a 
-b a 



After some algebra we have 

_a 2n 



(iph • J j 



. , , x 2 / _a 2n+l 

^(-ir +1 (n-j) , (Vj) 



(A.9) 



ip 2n+l (-l) n (ipn ■ j) , (A. 10) 
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and we may write 



exp (^pn ■ J j = 1 + (1 — cosy?) (h ■ + sirup (n ■ . (A. 11) 



The calculation of ■ J j is given by: 



(n-j) 2 = 



a 2 — 1 ab ac 
ab b 2 -l be 
ac be c 2 — 1 



Substituting equations (A. 9) and (A. 12) in (A. 11) we get 



(A.12) 



exp ( tph ■ J 



(1 — a 2 ) cos ip + a 2 ab (1 — cosy?) — c sin tp ac (1 — cos tp) + b sin tp 

ab (1 — cos ip) + csin y? (1 — b 2 ) cos <p + b 2 be (1 — cos ip) — a sin ip 

ac (1 — cos y?) — 5 sin y? 6c (1 — cos tp) + asmtp (1 — c 2 ) cos y? + c 2 

(A.13) 

Comparing (A.13) with the product 



R 3 (0)Ri(-0) 



cos — sin cos 9 — sin sin 9 
sin cos cos # cos sin 9 
— sin 9 cos 6* 



we get the following expressions for siny?, cosy? and components of ft: 



/ 20 2# 

sin 09 = 2 cos — cos —\ sin — 1 
T 2 2V 2 

COS + cos + COS cos 0—1 

cosy? = . 



• 2^-20 

sm" sin - sin — , 

2 2 2' 



(A. 14) 



(A.15a) 
(A. 15b) 



sin 9 (cos + 1) 



b = 



A <t> 9 l~- 2~d> , ■ 2~0 ■ 2~0 ■ 2 d> ' 

4 cos | cos | W sm | + sm | — sm | sm | 
sin sin 

4 cos | cos § J sin 2 § + sin 2 § - sin 2 f sin 2 \ 



sin (cos + 1) 



/I t) I ■ 2 4> i - 26* 2 6* • 2 tA 

4 cos | cos | W sin | + sin | — sm | sin | 



(A.16a) 
(A. 16b) 
(A.16c) 
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This result allows us to calculate [R 3 (0) R x (— 9)] N . 

r / \ t N 
[R 3 (0)Ri(-#)f = [exp (iph-3) 



cxp 



(Ntpn ■ j) . 



(A.17) 



The exponential exp will have the form identical to the matrix (A. 13), but 

with the substitution of ip by Nip. The components a, b and c, as well as simp and cosy? are 
shown in (A. 15) and (A. 16). 
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